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Abstract 

The phase transitions in Bose gases at constant volume and constant pressure 
are considered. New results for the chemical potential, the effective Landau- 
Ginzburg free energy and the equation of state of the Bose-Einstein condensate 
in ideal Bose gases with a general form of the energy spectrum are presented. 
q . Unresolved problems are discussed. 
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The Bose-Einstein condensation (BEC) is a phenomenon due to the quantum statistical 
correlations (pseudo-interactions) in ideal Bose gases (IBG) of non-interacting bosons. 
This is a condensation in the momentum space (Hk) but also it possesses some features 
of the usual (Van der Waals) condensation PJ El Ej This phenomenon strongly 
tJ- ■ depends on the spatial dimensionality d, and the energy spectrum e(k) = (h 2 k a /2m) 
of the bosons [0 < a < 2; k = \k\; k = {kj = (2n / nj / Lj)} , where j = l,...,d; 
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rij = 0, ±1, ...] in IBG in volume V = (Li...Ld) ~ L d and periodic boundary conditions 
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The original Bose-Einstein condensate (BEC) is not superfmid but the latter state is 
possible in many-body systems of interacting bosons (nonideal Bose gas, or, shortly, 
NBG) Pfl El El IH E] • The phenomena of BEC and superfluidity have a number of similar 
features and both of them are widely discussed in various problems of astrophysics 
(see, e.g., Refs. [3 III El IS] • In particular, BEC and superfluidity are relevant to the 
treatment of the so-called plasma-solid transition in "astrophysical matter" p], and 
for the behaviour of the neutron component of the dense matter in the interior of 
neutron stars 01313111; for applications to cosmological models of dark energy and 
dark matter, see Ref. [TUj . 

The main properties of BEC are known (see, e.g., Ref. |H ED U21 UH1 UH US] ■ Recently, 
BEC at constant density and constant pressure has been reviewed in Ref. ^EJ- Here we 
shall discuss the equation of state of IBG and BEC for spinless bosons. The treatment 
can be generalized for bosons with spin fT7\ IT%] . 

The superfluidity and the effect of interparticle interactions on BEC can be treated by 
both the mean-field like Gross-Pitaevskii ^J] and the microscopic Beliaev-Popov [2*Uj 
approaches. We shall discuss the latter within the renormalization group theory [3] 



I 



and for this reason we shall consider the following action of NBG (31 IH 120] : 

5[0] = -^G o 1 (g)|^(g)| 2 -^ £ P (qiWfaM<hMqi + to ~ ffs) , (1) 

<? 9i, 92,93 

— # 

where g = (a;/, fc) is a (d + l)-dimensional frequency-momentum vector, tu/ = 2nlkBT/h 
is the (Bose-)Matsubara frequency (I = 0, ±1,...), if)(q) is a C-number Bose field, u 
is the interaction constant, iV is the number of particle and (3 = 1/ksT. The (bare) 
correlation (Green) function Go(q) = (\if)(q) | 2 )o is given by 

G \q) = iut + eW+r, (2) 

where (— r) = /x < is the chemical potential of IBG (u = 0). 

Note, that the thermodynamics of NBG is given by the grand canonical thermodynamic 
potential 0(T, V, /x) = — /3 _1 Z, where the grand canonical partition function i£(T, V, /x) 
is defined by Z = J £>0exp (S) - a functional integral over the possible field configura- 
tions. The usual field theoretical investigations of the action (1) are performed in the 
thermodynamic limit: iV — > oo, V — > oo, provided < p = (N/V) < oo. 

Using the models of IBG {y = 0) and NBG defined by Eqs. (1) - (2) we shall consider 
the equation of state for BEG For IBG this can be performed exactly, whereas for 
NBG we must use the loop expansion (see, e.g., Ref. [3 ). The main features of BEC 
can be revealed within the one-loop approximation [SI EI]- We shall discuss the effect 
of the inter-particle interaction on BEC by using results from preceding works. 

We wish to emphasize that the phase transition to BEC cannot be easily put to the 
usual classification of phase transitions 3\. Above T c this phase transition resembles 
certain features of second order phase transitions but is rather different from the stan- 
dard notion about these transitions 0. On the other side, the equation of state of 
IBG below T c is quite similar to known equations of (almost) first-order phase transi- 
tions and tricritical points |3] ■ There is a close similarity between the phase transition 
properties of IBG and the spherical model in the ferromagnetism, and this point will 
be discussed in the remainder of this report. In our consideration we shall essentially 
use results from preceding works jTHl El El El El E] (for a comprehensive reference 
to original papers, see, e.g., the reviews HI)- We work with general values of d and 
a but our consideration includes the important case of d = 3, a = 2. 

2. Free energy and thermodynamics 

For a convenience we shall introduce a fictitious external field h which is thermo- 
dynamically conjugated to the order parameter \1/ = ("0(0)) of uniform BEC. Note, 
that in case of uniform BEC the (q = 0)-mode ^(O) can be represented by the sum 
^(0) = (^(O)) + Sip(0), where 5i[)(0) is the (uniform) fluctuation mode. The consider- 
ation of an external uniform field h can be performed by adding a term 

S h = i N>*(0) + c.c.)} . (3) 
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The thermodynamic potential of IBG can be written in the form 

fi(T, r, h) = -r'VX^Aid, a)g d/a+1 (Pr) - (4) 
where g v {y) is the Bose function [3JE3EI|> the thermal wavelength is given by 

and 

The potential (4) obeys the differential relation dVt = -SdT + Ndr - ^dh* - ^*dh. 

By a suitable Legendre transformation we obtain another thermodynamic potential 
Q(T,r, where the natural variable is 

il = -r'VX^Aid, a)g d/a+1 ((3r) + Nr^ 2 . (7) 

We can restrict ourselves, without a loss of generality, to real h and \&. The susceptibil- 
ity is then xt — d^/dh = 1/Nr ~ t -7 , where 7 is the susceptibility critical exponent 
and t — (T — T c )/T c (see,e.g., Ref. 0). 

For small r and energies (e ksT), which correspond to the critical regime near 
the phase transition point, the correlation function x{k) = Gq(0, k) takes the form 
~ (ck a + r). This gives us the correlation length £ = (c/r) 1 / " ~ t~ u and 
x(k) ~ k~ 2+v at r = 0. From x{k) ~ k u we obtain that the Fisher exponent 77, defined 
by x{k) ~ k~ 2+v , is equal to (2 — a) for all dimensional ranges and possible constraints 
(of constant volume V or constant pressure P). 

In order to obtain the correlation length exponent v and the exponent 7 of the suscep- 
tibility xt, we need to calculate the function r(t). The latter is different for the cases 
V = const and P = const. In the thermodynamic limit, the constraint of constant 
volume is equivalent to a constraint of constant density (p = const). The results for the 
critical exponents corresponding to these cases are summarized in Ref. ^B] for various 
spatial dimensions d. Here we restrict our consideration of the critical exponents and 
the equation of state to features which demonstrate the difference between the phase 
transition to BEC and both first and second order phase transitions for the most in- 
teresting interval of spatial dimensions a < d < 2a that includes the case a — 2, d — 3. 



3. Constant density 

To consider the effect of the thermodynamic condition of constant density p = (N/V) 
we need the free energy F = Vf(T,p,^f). To obtain this thermodynamic potential 
near the phase transition point to BEC we expand the Bose function in (4) in powers 
of j3r and express r as a function of p. This procedure is performed for the potential 
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Q given by Eq. (7). Further we obtain the potential F with the help of a Legendre 
transformation of the form: 

F = Cl- Nr\ r=r{p) , (8) 
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on 

Vdr 



X T d A{d ) <j)g d/a {{3r) + p^ 2 . 



(9) 



The lowest temperature T c , at which the system does not condensate (\I/ = 0) is 
obtained from Eq. (9). When d < a, the Bose function is divergent for (/3r) — > which 
means that BEC may occur only at the absolute zero (T c = 0) - a zero temparature 
BEC PSJ. When d > a, we obtain that HEU 
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> o, 



(10) 



where g d / a (fy = ((d/a) is the zeta function. Below the critical temperature T c , BEC 
occurs (^ > 0). 

For the most interesting case of dimensions a < d < 2a, the result for the free energy 
density f(T,p, \1>) to the lowest order in \t\ -C 1 is given by: 



where 



C f 



d 

- - 1 

a 



C(d/a) 



d/(d-a) 



a/(d-a) 



(k B T c )p . 



(12) 



\T(l-d/a)\_ 

In our derivation of the energy (11) a ^/-independent term has been neglected. Such 
terms can be ignored because they belong to the energy of the disordered phase (\& = 0). 
For this reason the net free energy of the BEC can be obtained from (11) by neglecting 
the ^/-independent term of type t d ^ d ^ a \ The free energy (11) is of Landau- Ginzburg 
type 0; one may easily expand /(^l/) in powers of \I/. 

The external field h is fictitious and can be neglected in studies of the thermodynamics. 
In this case, the equation of state h ~ d(f /d^) becomes (df/d^) = 0. From Eq. (11) 
we easily obtain the equation of state for h = 0: 



* ( ^ 2 + - t 

a 



cr/(d-o) 



0. 



(13) 



The solutions of Eq. (13) are: \1/ = (disordered phase), and \1/ 2 = —{d/a)t > 
(corresponding to BEC for t < 0). Note, that the chemical potential (p = —r) can be 
obtained in the form 
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-1 0-/(d-cr) , 

,,,, , u (* a + V/C*-). 

|r(l — a/0") | J a 



ad/a) 



(14) 
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BEC is possible for t < under the condition /i = 0. The latter introduces thermody- 
namically forbidden (unstable) domains in the phase diagram. 

The main conclusion which can be drawn from the Ginzburg-Landau free energy (11) 
of IBG at constant p is that the point t = resembles a trictitical point 3J. Usually 
such multicritical points occur on a phase transition line where the phase transition 
changes from second order to a symmetry conserving first order phase transition (or 
vice versa). Here this is not the case but the similarity with the usual tricriticality is 
in the fact that the coefficients of both the \1/ 2 — and the \1/ 4 — terms in (11) tend to 
zero for t — > 0. On the other side the phase transition to BEC is a continuous phase 
transition which exhibits critical exponents identical to the critical exponents known 
from the Berlin-Kac spherical model in ferromagnetism |22| 128] (see also Refs. 0I3])- 
Here the condition of constant density p plays the role of the spherical condition for 
the spins in the spherical model. 

Interaction effect. Most of these exceptional properties of the phase transition to BEC 
in IBG are not present in systems of interacting bosons described by NBG. In this case 
an additional term of type w\l/ 4 will appear in the effective free energy f(T, p, and the 
respective coefficient of the \l> 4 — term will remain finite at T\ - the critical temperature 
to superfluid state in interacting Bose fluids (gases and liquids). The generalization of 
our treatment to the case of interacting bosons described by the action (1) may lead to 
new thermodynamic properties. A similar generalization for the spherical model was 
made in Ref. Another way of treatment of NBG may be performed within the 

loop expansion [21]. It is believed that the A— transition described by NBG belongs to 
the so-called XY universality class of standard second order phase transitions 

Recent studies of interacting bosons indicated discrepancies in the theoretically pre- 
dicted values of the critical temperature T\ (231 |2Hl EH1 EZj . The problem for the 
calculation of the phase transition temperature of interacting many-body systems is 
a hard and still unresolved problem of the theory The fluctuation shifts of the 
critical temperature that are usually calculated from field models (see, e.g., Ref. [3| are 
very small and do not include essential contributions due to the large-momentum (high- 
energy) fluctuations of the order parameter field ip. Therefore, the correct treatment of 
the phase transition temperature in many-body systems with interparticle interactions 
requires new theoretical methods. 

4. Constant pressure 

The only papers where the effect of the constraint of constant pressure on the phase 
transition properties of IBG has been investigated so far are Refs. [T51 IT41 H[ I16j. This 
case should be taken in mind in interpretations of real experiments, in particular, in 
low-temperature experiments on BEC in trapped atomic gases (see. e.g., Ref. [TH]). 
where the density p varies but the pressure P is (almost) fixed. In this case BEC 
occurs at finite temperatures (T c > T > 0) for all spatial dimensions d > 0. The 
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critical temperature will be jl] 



T C (P) 



v/{d+o) 



(15) 



_C(1 + d/a)A(d,a)k B _ 

where Ao is given by A = AT 1 /' 7 and Eq. (5). For d = a = 2 one obtains the result for 
T c known from Ref. [TBI . 

Although a number of results are known from preceding works, in particular, for the 
two-dimensional case (d = 2) the entire picture of the phase transition to BEC at 
constant pressure is not still clear for all spatial dimensions d. In particular, this is the 
case of interacting bosons. We have no information about research papers devoted to 
this problem. 

5. Concluding remarks 

We have presented a brief discussion of several properties of the phase transition to 
BEC. New results have been obtained for the effective Landau-Ginzburg free energy 
(11), the equation of state (13) and the chemical potential (14) of BEC in ING for 
general values of d and a. 

The effects of the constraints of constant density and constant pressure on BEC in ING 
and the condensation to a superffuid phase in systems of interacting bosons (NBG) are 
not yet clarified in a comprehensive way. The renormalization group methods j3J B] 
reveal a dimensional (quantum to classical crossover) in interacting Bose systems and 
other basic universality features of the so-called quantum phase transitions at zero and 
very low temperatures but the equation of state below the phase transition point is 
not investigated in details. This is the situation for the whole variety of Bose systems 
known in condensed matter physics j3J E] and, in particular, for Bose fluids of real 
atoms. 

For Bose fluids the concept of universality of the quantum critical phenomena seems 
to be invalid and mainly for this reason the quantum phase transitions of second order 
can be classified in two groups: as universal and non-universal (see Ref. jl]). However, 
a number of examples indicate that the quantum phase transitions are often of first 
order or are described by multicritical points which are different from the critical points 
of standard second order phase transitions. 

In astrophysics, a quantum phase transition (T c ~ 0) may occur in cases of very low 
density or very low pressure of the respective Bose fluid; see Eqs. (10) and (14). For 
the matter in the interior of neutron stars one should investigate the so-called classical 
limit jl] in which the quantum fluctuations are irrelevant. This is the case of T c > 
discussed in the prevailing part of our report. 

BEC of spin bosons ^7] can be treated within the framework of a generalization 
of our treatment. For this aim one may consider a complex vector field ip(q) = 
{ipal « = 1, ...n/2}. For n = 2 one obtains the complex scalar field in Eq. (1). When 
such a system is placed in an external magnetic field, essentially new results can be ob- 
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tained fT7j . The unresolved problem about the superfluidity of interacting spin bosons 
in external magnetic field is of special interest in astrophysics. 
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